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Abstract 



We analyze the fully relativistic, field-theoretical treatment of the scalar Coulomb problem. We work 
in a truncated Hilbert-Fock space containing the two-constituent states and the two-constituent-and-one- 
massless-exchange-particle states. Self-energy contributions are dominant for large values of the coupling 
^s^j ■ constant. Using the covariant formulation of the self-energy does not alter these results significantly. Both 

' the weak-coupling limit and the heavy-mass limit lead to the non-relativistic results. 



1 Introduction 

In this paper we will solve the simplest possible scalar field-theoretical Coulomb problem in a Hamiltonian 
framework B. This approach is much simpler than most other field-theoretical approaches, and yields better 
results. Furthermore, we also include the self-energy parts, which are usually ignored, and find that they 
give significant contributions. Our main motivation for this work lies in improving our understanding of the 
field-theoretical bound-state problem before we turn our attention to stronger binding, and more physical, cases. 



Qh! There is a wealth of relativistic bound-state equations |^|, ||, Q| . Most of them are three-dimensional reductions of 

the Bcthe-Salpeter equation 0] in the ladder approximation. Since there is a degree of arbitrariness involved in 
this reduction R |?j , the equations can be tailored to satisfy unitarity (|] , gauge invariance || , or the one-body 
limit ||, which are some of the typical problems with the underlying Bethe-Salpeter equation |)| |icj] . However, 
these are approximations upon an equation with problems, which is a questionable practice. 
Generally, it is difficult to favour one equation above the other jy] . Therefore, there has been some attempt to 
establish standard results [fijl f?J, which include all orders in perturbation theory, in the quenched approxima- 
tion. To do so one starts from the path integral for the two-particles four-point function and integrate out some 
degrees of freedom analytically, and the resulting integral is integrated numerically in the Feynman-Schwinger 
representation by the Monte Carlo method. The typical problems with this method are related to the discretisa- 
tion of space-time, the Wick rotation, which leads to decaying states of which only the lowest-energy states can 
be determined accurately, and the problem to handle massless particles satisfactorily [Q. However, the results 
are promising and show the underbinding of the Bethe-Salpeter equation in the ladder approximation |l3|, Q . 
Another, distinct approach is based on Haag's expansion of asymptotic states in field theory in free and bound 
states, which lead to an hierarchy of coupled equations in the coupling constant expansion (T5J . With some 
difficulty the scalar field theory, which is also analyzed here, is solved. The results seem to indicate a failure to 
obtain both the one-body limit and the weak-coupling limit in the Haag-expansion approach. 
Most prominently has been the return to Hamiltonian field theory (T| . Because the bound-state is a stationary 
state, and manifest covariance is lost, it seems appropriate to solve the bound-state problem in the Hamiltonian 
framework. Most calculations are done in the light-front Hamiltonian framework |n| [lTj] , which preserves part of 
the boost symmetry at the cost of rotational symmetry . We will use the ordinary instant-form Hamiltonian 
instead to show that that indeed yields the same results as the Kepler problem |is|| , in the weak-coupling limit. 
The reason to study the case of massless exchange particles in particular are manifold. Firstly, the massless 
particles are a typical and hard problem for field theory, as they generate long-range interactions, and lead to 
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singularities in the S-matrix associated with multiple soft-photon emissions [|o|, HJ] . None of these problems 
have to be present in the bound-state calculation. Indeed, in the Hamiltonian approach the infrared singularity 
is absent, contrary to the nonrelativistic treatment using the Coulomb potential and the covariant approach. 
When constituents are bound, they are off-shell, which means they do not share enough energy among them to 
become on-shell. Therefore they do certainly not have enough energy to radiate off soft photons, which is the 
origin of the infrared complications in the covariant approach J2^, |2^] . 

Consider the S'-matrix element of two constituents exchanging a massless particle. This matrix element has a 
singularity. It occurs because at zero momentum exchange the difference in energy of the initial and intermediate 
states vanishes. If the two constituents belong to a bound state, the singularity is hidden by the binding energy. 
Apparently, the binding energy has the same effect as the mass of the exchanged particle: screening at long 
distances. 

Secondly, eventually, one has to admit that all fundamental interactions that generate bound states are massless, 
gauge theories. For the case of positronium, a relativistic approach can be used (23), however, rather in the 
cumbersome Coulomb gauge |^4| where the lowest-order kernel reduces to an instantaneous interaction, and 
many other cancellations occur. For the heavy-light case, like the hydrogen atom, the starting point has to be the 
non-relativistic Coulomb potential J25| , or the equivalent Dirac equation with such a potential. This Coulomb 
potential has survived many revolutions in physics, mainly because of its phenomenological success: it predicts 



atomic spectra to a high accuracy |26], 25 1. This success stood in the way of fundamental development, where 
the classical action-on-a-distance concept is replaced by the correct notion of particle exchange . There 

is no viable alternative which reproduces the particular long-range interaction generated by the exchange of 
massless particles. At least, there is no alternative that does not introduce additional problems in the one-body 
limit. 

In this paper we recover the typical scaling associated with the weak-coupling Coulomb ground state. This 
paper is not just a formal exercise to link theory with experiment. Eventually, we like to understand how 
relativistic effects come about in deeply bound systems, and give sensible answers to where the spin of a system 
is coming from, how the radius of the bound state is related to the masses and the coupling strength, what 
the mass of a constituent is, and which are the effective degrees of freedom. We emphasize deeply bound states, 
i.e., states with a binding energy of the same order as the masses of the constituents, because we feel that for 
them the problems we want to address are most pressing. For weakly bound states one can treat relativistic 
aspects, like relativistic kinematics and the creation and annihilation of particles as a correction to the classical 
approach, where a potential picture is appropriate. Although it would be nice to understand better the origin 
and limitations of the potential picture. 

The Hamiltonian framework is conceptually quite different from the standard, covariant field-theoretical frame- 
work, although the final formulae may be similar. In this framework exchange diagrams exist, but not as a 
potential. Rather, they are to be interpreted as matrix elements of the fundamental, local interactions, between 
different Fock sectors. 

In the fundamental description of bound states there is no place for the concept of a potential. All interactions 
are local. What seems to be a non-local interaction is due to the smearing of the constituent particles, which 
are not just single particles but a linear combination of states with different numbers of bare particles. 
The Hamiltonian description gives us a consistent way to truncate the infinite complexity arising in field theory, 
by truncating Fock space to two- and three-particle states. Many arguments concerning bound-states are 
devoted to the truncation to a finite number of particles, since, in order to generate a bound state which 
arises only non-perturbatively, an infinite summation of diagrams is required. So the question is which finite 
set one should take to iterate infinitely many times. For the Bethe-Salpeter equation |j| the results depend 
crucially on the answer to that question. However, formulated as a Hamiltonian eigenvalue problem this infinite 
summation appears naturally, given the Fock-state truncation. If this truncation to the two- and three-particle 
Fock-state sector is legitimate anywhere, it should be here, because in the weak-coupling limit a few particles 
are expected to contribute and anti-particles (Z-diagrams, which appear in a covariant approach) are expected 
to be suppressed, due to their large virtuality. We shall analyze the validity of these arguments. 
So, for our problem we assume that our bound state consists of two pieces: one with two oppositely charged 
scalar particles and one with an additional, massless scalar particle that can be emitted and absorbed by each of 
the charged particles. We use scalar particles to keep the problem as simple as possible, and free from questions 
associated with gauge symmetry. A careful analysis of the three-particle state is required in order to recover all 
of the physics. One finds an exchange contribution and a self-energy contribution; the latter is infinite and has 
to be renormalized M. Of course, we face the problem of the interpretation of divergent integrals. 
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The concept of self-energy is still a mystery to many. In solid-state physics one can accept that a moving 
electron deforms the crystal, generating a drag which can be interpreted as an effective mass. In deep-inelastic 
scattering one can interpret the self-energy as the mixing-in of multi-particle states as sufficient energy becomes 
available when the particle goes off-shell. In a bound state the missing energy is the binding energy, which does 
not belong to a particular particle, but affects each of the particles separately. One may look upon this in the 
following way: each constituent is not a single particle; partly it is one particle, partly it is two, and so on. 
The relative weight of each of these states of a constituent is determined by two things: the number of states 
available and the virtuality of these states cjj. In a classical world only the one-particle state wq would exists 
because that state has the lowest energy luq < u>i. In the quantum world and within the Hamiltonian framework, 
energy is only conserved asymptotically, so higher energy states, when available, will also be occupied. The 
weights of virtual states is determined by their virtualities. If we now put a "compound" constituent in a bound 
state with binding energy Eb, the relative virtuality u3i/u>o > (cJj + Eb)/(u>o + Eh) of the multi-particle states 
in the constituent is less than for the free particle, therefore the summed energy will lower, and the self-energy 
will increase the binding. 



2 Theory 

In the present paper we solve the scalar Coulomb problem. We limit ourselves to this case in order to avoid 
numerous problems concerning gauge theories. Here we are rather interested in the conceptual problem of what 
a bound state is and what determines its properties in the simplest possible case. Spin effects and questions 
concerning gauge invariance can be dealt with later. 

Our model Lagrangian describes the interaction of three scalar fields <p a , 4>b, and X , 

£ = ~i>an + ZaX]4>a - \^d^ + m\ + e bX ]^ b - \x^X , (1) 

which have the respective masses: m a , nib, and zero. It will allow for a bound-state solution, if the charges, 
e a and eb, with which the massive fields couple to the massless field x are imaginary and have opposite signs. 
This choice makes the Lagrangian non-Hermitian. However, also with real coupling constants this Lagrangian 
would be "sick," since then the spectrum is not bounded from below pjl . However, for a bound state we 
restrict ourselves to low orders in particle number, or coupling constant, which, we will show, leads to sensible 
results, apparently unaffected by the formal problems with the Lagrangian. This Lagrangian would lead also 



to the Wick-Cutkosky model |30| as only the product e a &b would appear. The only difference is the opposite 
sign of the self-energy part, which is not present in the original W-C model. For real couplings of the same 
sign, the problem is similar to scalar gravitation where the same charges attract, which is precisely the reason 
why the vacuum is unstable, as more and more charges together have a lower and lower energy. The opposite 
imaginary charges would correspond to electrodynamics, which would, if implemented properly, require vector 
coupling and gauge fields. However, to recover the Balmer spectrum this is not necessary, and we want to keep 
the discussion clear by avoiding issues concerning gauge symmetry. However, a more important point to make 
is that, unlike in the W-C model, the heavy-mass limit, or one-body limit, where one constituent has a infinite 
mass yields the correct solution [^7], [2^]. Although the spectrum of the W-C model yields the Balmer series 
in the weak-coupling limit, the equation does not reduce to an effective one-body equation, where the heavy 
particle plays the role of a static source. 

The transition from the Lagrangian formalism to the Hamiltonian one is rather tedious, but can be found in 
the literature [^7], EL Here we only like to make some general comments. In order to give the theory a particle 
interpretation the free fields are quantized. Therefore, if a field f(k) is known on a space-like hypersurface in 
four-dimensional space-time, it is known in the whole space, via the mass relation: 

-f— /(w ro (k),k), (2) 

where the phase-space factor (f> m depends on the mass m and the three-momentum k perpendicular to the 
time direction r. The on-shell energy w m (k) is the solution of the equation k 2 — m 2 = 0. The 8 T function 
restricts the integral to the positive energy solutions. Negative energy solutions are associated with holes or 
backward-in-time moving particles, which, in the Hamiltonian formulation, are the charge- and parity-reversed 
partners, and should be treated as distinct particles. 
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Each intermediate particle is associated with a phase-space factor which, for a time-reversal invariant 

formulation, should be divided among the creation and annihilation vertices of this particle. We will work in 
the ordinary equal-time formulation, so r = x° and 

w m (k) = Vm 2 + k 2 , (3) 
4> m {k) - 2w m (k). (4) 

For the simplest treatment of the field-theoretical problem above we assume that the bound state consists of 
only two different Fock states: 

\i>) = \PaPb) + |k a k fc k 7 ) . (5) 

No further approximations are required. The Hamiltonian restricted to this subspace of the Hilbert space allows 
only for the possible absorption or emission of one massless particle: 

where the diagonal entries are the on-shell energies of the respective Fock states, and the off-diagonal entries 
are the absorption and emission of the \ field: 

w 2 = u a (p a ) + w b (p b ) , 

< = 

u 3 = w a (k a ) + w fc (k b ) +w 7 (k 7 ) . (7) 
Here we introduce [• • •] as the shorthand notation for the normalized Hamiltonian interaction 

rk a k b -| = f d 3 p c S 3 (k a +k b -p c ) ^ 

r P . 1 f d 3 k c d 3 k d 5 3 (p a - k c - k d ) 

LkckdJ " J ^4> m APa)&nM&nM ' 

and {• • •} generates the proper change of variables for the particles not affected by the interaction: 

{ll}= Jd 3 p b 6 3 (k a -p b ). (10) 

We solve the eigenvalue equation above by introducing the eigenvalue £asa Lagrange multiplier, and expressing 
the three-particle state in terms of the two-particle state and a function of the unknown E, from the eigenvalue 
equation: 

lkkk , <:}fc k i+<}[p: k i fu , 

k a k 6 k 7 ) = — 7^ 77—^ ~— r PaPb) • (11) 

E - w a (k a ) - Lo b (k b ) - w 7 (k 7 ) 

This three-particle state shows the interesting features of a bound-state problem. Although a third particle is 
present, and so the number of degrees of freedom has increased, the actual degrees of freedom are still those 
of the two-particle wave function, and an additional parameter in the form of E. This wave function has 
significant values only in restricted areas of the k a ® k b ® k 7 space; however, large back-to-back values can 
occur, unquenched by the two-particle wave function. 

The dimensionful coupling constants ie a / \/m a m b = —ie b / 'y/m a m b = e are related to the fine-structure constant 
a = e 2 /(167r), in such a way that the bound state is neutral. In order to find the bound state we use two sets 
of variational wave functions, tuned by the parameter /j,: 

|P, -P>cov : 



^/Lu a Lu b {p - m a - m b + uj a + Lu b ) 2 

1 



|P '- P> - = (2m M + P 2 ) 2 ' (12) 
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Figure 1: A graphical representation of Eq. |IJ 



where m is the reduced mass. We will work in the center-of-mass frame of the two constituents, therefore, 
uj a = Lu a (p) and uj b = uj b (—p) = w&(p). The second wave function is motivated by the non-relativistic ground- 
state wave function, while the first wave function would follow from a relativistic vertex in Hamiltonian form, 
after it is differentiated with respect to the scale parameter ji. In practice the form of wave function does not 
matter, although the actual value of the parameter /i does. The energies differ only 0.1% for the different choices 
of wave functions. 

For each bound-state energy E the coupling constant a is determined by minimizing it with respect to the 
wave-function scale parameter jx: 



E(p, -p|p, -p) - (p, -pMp, -p) 
a _1 (p, -p| < |k a k&k 7 ) 



(13) 



Note that the quantity in the denominator is actually independent of a. It contains two types of contributions: 
exchange terms and self-energy terms, 



(p, -p| < |k a k b k 7 



(V ab (E)) + (V ba (E)) - ( 



(14) 



where (• • •) indicates the expectation value for the wave functions ip of Eq. |12|. Equation 14 is graphically 
represented in Fig. |l|. The exchange part is finite 



(Vab(E)) = 



am a mb 



d 3 k 



Hp - k) 



4tt 2 J y/u a (p)u a (p) J ^/uj a {p - k)uj b (p - k)w 7 (k) E - uj b {p) - w a (p - k) - w 7 (k) 



however, the self-energy part is divergent: 



£a 



am a m b 



d 3 k 



An 2 J uj a (p - k)w 7 (k) E - uj b (p) 



(p-k) -w 7 (k) 



(15) 



(16) 



The integral is divergent because of the occurrence of the forward part of the one-loop self-energy correc- 
tion. We renormalize this contribution by subtracting the on-shell value: S r (_E) = 'S(E) — E(w a (p) + uj b (p)). 
Rcnormalization means the redefinition of the parameters, like masses and charges, of the theory. 
We use the on-shell renormalization J2£| , which boils down to discarding any change in the values of the physical 
quantities for the case of free, constituent particles. This will give the same results in the same approximation, 
for observables, as any other approach, but frees us from lengthy discussions on counterterms and such. However, 
there will be small differences in the results, depending whether on-shell means on the mass-shell [j? = m 2 ), 
as in the covariant approach, or on the energy-shell (E — uj), as advocated here in the Fock-state truncation. 
The angular integrations can be performed analytically. In order to calculate the integrals with the least 
numerical problems we integrate the self-energy contributions analytically, thereby avoiding the integrable 
singularity in the subtraction term S(w a + u> b ) in our numerical calculations. For the sake of completeness we 
give the result of a tedious but straightforward calculation: 



fcad/c a 



In 



LU b - E + LU a + \k b + fcqj 

LU b - E + LU a + \k b - k a \ 



m 2 a ln[m a ] 



subtracted 

u b — E 



1 



1 



U) b - E k b — LU a 



Ua 

ml (E - uj b + LO a )(E - uj a - Lu b ) + ml 

™2 



ml ln[fc b + uj a 

\ In[w„ + LJ b - E]{u a + Ub -E)(l + ": a 

2 \ (pj b - E - k b ){k b + Lo a ) 

- In /,'/, -!- j.; h - E ■ 



{2k b + u a + u b - E) 1 + 



(h + uj b ~- E)(k b + Lo a ) 



(17) 
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Figure 2: The properties of the ground state compared to the non-relativistic results, in the case of the masses 
relevant for the hydrogen atom, m a = rriy^m,^ = m e . 



The integrals converge well. With a typical compact integration coordinate the integrands are smooth bell 
curves in the middle of the domain, from which we know that we do not miss any contributions for small or 
large values of the momentum. We use a simple integration procedure and with 500 integration points we 
achieve a sure accuracy of 4 to 5 decimal places, which is in the same region as the systematic error due to 
the wave- function Ansatz. It is interesting to observe how different the self-energy contributions are for the 
heavy-mass and the equal-mass cases, shown in Figs. ^ and [| respectively. 

Since Hamiltonian renormalization is not a well-established procedure, although we are convinced of its legit- 
imacy, we repeated the same calculation with the covariant expression for the self-energy, which means that 
only for that part we also include the states with three particles, one anti-particle, and a photon, which leads 
to the known covariant result for the self-energy: 



o 

TO 



2 

P 



(18) 



where £ c (p 2 = to 2 ) = 0, which corresponds to the on-mass-shell renormalization; for the physical mass there 
is no mass correction. For small coupling the result differs only little from the restricted Fock-state result with 
m a = nib = to. For larger values of the coupling constant the self-energy contribution is slightly suppressed in 
the covariant approach, as can be seen in Fig. 

The heavy-mass limit yields the correct result, which reduces to the non-relativistic result in the weak-coupling 
limit. The latter problem is dealt with in a rather cumbersome way in the Bethe-Salpeter equation and related 
approaches p7|, as the one-particle exchange approximation (the "ladder approximation") does not work there. 



3 Renormalization 

Since we solve only the lowest non-trivial case of binding with a massless exchange field, there are many 
contributions we do not include. Eventually, to show that the lowest order includes most of the relevant physics 
of a bound state, we have to pursue the problem of convergence. This is part of the future work. However, 
since renormalization involves divergent contributions already present at the lowest order, also here some care is 
required. In the renormalization different orders in the coupling constant are mixed, while in the Hamiltonian 
approach some of these terms, which should be grouped together, might not appear. 

For example, from the Ward identity, it would follow that one should combine wave-function renormalization 
of the self-energy diagrams with the vertex correction, to yield a finite answer. The former diagram appears 
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Figure 3: The properties of the ground state compared to the non-relativistic results, in the case of two equal 
masses, m a = rrib = m. 




Figure 4: The comparison in the equal-mass ground-state binding energies and the self-energy contributions 
to that result between the case where the self-energy is renormalized on energy-shell, which includes only the 
forward diagram, and the case where one-particle mass-shell renormalization, from the corresponding Feynman 
diagram which also includes the anti-particle contributions, is used. 
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in our equation, but the latter do not. We choose not to perform the wave-function renormalization, for one 
important reason: it leads to a trivial theory, with zero coupling, in our case of massless exchange particles. 
In the case of finite exchange-particle masses the renormalized coupling constant scales with the logarithm of 
the mass. Only if the vertex correction is also included, the infinite wave-function renormalization is counter- 
acted by the infinite charge renormalization, and only a mild dependence is left. Since these divergences are 
infrared divergences, it is even more physical |p2|| to group the wave- function renormalization with the vertex 
renormalization, and have both or none. In our case it is none. 

If we compare our result with the work by Ji |32| we find, as expected, the opposite effect from the self-energy, 
because we have chosen the starting Lagrangian in such a way that we have opposite self-energy contributions; 
in our Lagrangian identical particles repel each other, like they should do in the Coulomb problem. However, 
unexpectedly, one can conclude from Ji's work that the self-energy cures a bit of the illness of the <fi 3 theory. 
One would expect the binding energy to make a nose dive as the coupling constant increases, due to the 
unboundedness of the Hamiltonian. However, including the self-energy decreases the binding in the case of the 
scalar exchange as studied by Ji. 

For the Bethe-Salpeter equation, with dressed propagators, it is found that the self-energy increases the binding 
slightly, however, these results are for a relatively large mass of the exchange particle, and the wave-function 
renormalization is used, which counteracts the weaker binding caused by the mass renormalization[Q. More- 
over, other problems appear in the covariant Bethe-Salpeter formalism, such as unphysical solutions |33"[ |To| , 
which have, at least in part, to do with the relative time-coordinate appearing in the equation. However, there 
are claims that these problems result from the unphysical Lagrangian |34). For the moment we want to conclude 
that if renormalization is pursued beyond the simplest on-shell subtractions, many questions are still open. The 
same conclusion can be drawn from the work of Glazek et al.]35|] on light-front Hamiltonian Yukawa theory, 
which seems, insofar, to confirm Ji's findings for self-energy effects of real scalar exchange fields. However, it 
must be noted that in all these cases the Lagrangian was different from our Lagrangian, which in practice boiled 
down to the opposite sign for the self-energy. This is altogether a different problem. When we reversed the sign 
of the self-energy in our calculations, as a simple check, we found instabilities for a > 0.5. This is a subject for 
further study. 

4 Conclusion 

Generally, one can say that the approach advocated here is straightforward compared with the Bethe-Salpeter 
equation and subsequent three-dimensional reductions. The normalization is that of ordinary quantum mechan- 
ics, and spurious solutions, associated with excitations in the relative time do not occur, since the Hamiltonian 
framework has only one reference time. Moreover, issues concerning anti-particle content of the wave function, 
present in the covariant formulation, do not arise here, since the approach is based on the Fock-state truncation. 
Another advantage over the standard covariant approach is the absence of infrared singularities, due to the fact 
that the states are always off-energy-shell, and these singularities are associated with on-energy-shell states with 
vanishing photon momentum. 

Apart from the technical advantages of this approach it is important to notice that the self-energy contributions 
dominate the ground state in the strong binding regime. If we perform the same calculations including all the 
relativistic effects except for the self-energy contribution, we find completely opposite results | BjJ ; the binding 
energy is lower than the non-relativistic one, as can be seen in Fig. ^ This indicates that the usual description 
of strong-coupling binding in terms of exchange diagrams is, to say the least, incomplete. 

5 Acknowledgments 

We like to thank Prof. John Tjon, Prof. Reinhard Alkofer, and Prof. Franz Gross for remarks and discussion 
regarding this work. We also like to thank Jim Friar, Roland Rosenfelder, and Axel Weber for their comments. 

References 

[1] W. Heitler, The quantum theory of radiation, (Clarendon, Oxford, 1954). 
[2] I.T. Todorov, Phys. Rev. D 3 (1971), 2351. 



8 



-0.5 



LLJ 
UJ 



■1.0- 



// 

/ 

t 



— with self-energy 
-- without self-energy 
-- Wick-Cutkosky 

- Light-Front 




l i i i i i i i l 

0.0 0.5 1.0 1.5 2.0 

a 

Figure 5: Comparison between the result for the full calculation, and when the self-energy terms are dropped 
for the equal-mass case. We also compare the results with the Wick-Cutkosky model |3C| [l(J and light-front 
results |3l) . 

[3] R. Blankenbecler and R. Sugar, Phys. Rev. 142 (1966) 1051. 
[4] F. Gross, Phys. Rev. 186 (1969) 1448. 

[5] H.A. Bethe and E.E. Salpeter, Phys. Rev. 82 (1951) 309. E.E. Salpeter and H.A. Bethe, ibid. 84 (1951) 
1232. 

[6] F. Gross, Phys. Rev. C 26 (1982) 2203. 

[7] T. Nieuwenhuis, PhD thesis Utrecht, Oct. 1995. 

[8] D.R. Phillips and S.J. Wallace, Few Body Syst. 24 (1998) 175. 

[9] C. Itzykson and J.B. Zuber, Quantum Field Theory, (McGraw-Hill, New York, 1980). 
[10] N. Nakanishi, Prog. Phys. (Suppl.) 43 (1969) 1; ibid. 95 (1988) 1. 

[11] C. Savkli, J.A. Tjon, and F. Gross, Phys. Rev. C 60 (1999) 055210; Erratum-ibid. C 61 (2000) 069901. 
[12] Yu.A. Simonov and J.A. Tjon, Ann. Phys. (N.Y.) 228 (1993) 1. 
[13] T. Nieuwenhuis and J.A. Tjon, Phys. Rev. Lett. 77 (1996) 814. 

[14] C. Savkli, Proceedings of the European Conference on Problems in Few-Body Systems, Evora, 2000 (to be 
published). 

[15] OW. Greenberg, R. Ray, and F. Schlumpf, Phys. Lett. B 353 (1995) 284. 
[16] S.J. Brodsky, H.-C. Pauli, and S.S. Pinsky, Phys. Rept. 301 (1998) 299. 
[17] S.J. Brodsky, J.R. Hiller, and G. McCartor, Phys. Rev. D 60 (1999) 054506. 
[18] P.A.M. Dirac, Rev. Mod. Phys. 21 (1949) 392. 

[19] S. Fliigge, Practical quantum mechanics, (Springer, Heidelberg, 1971). 
[20] F. Bloch and A. Nordsieck, Phys. Rev. 52 (1937) 54. 



9 



[21] S. Weinberg, The quantum theory of fields, Vol. 1, (Cambridge, Cambridge, 1995). 

[22] J.M. Jauch and F. Rohrlich, The theory of photons and electrons, (Springer, New York, 1976). 

[23] T. Murota, Prog. Phys. (Suppl.) 95 (1988) 46. 

[24] G.S. Adkins, P.M. Mitrikov, and R.N. Fell, Phys. Rev. Lett. 78 (1997) 9. 

[25] M.I. Eides, H. Grotch, and V.A. Shelyuto, hep-ph/0002158, and the references therein. 

[26] H.A. Bcthc and E.E. Salpeter, Quantum mechanics of one- and two-electron atoms, (Springer, Berlin, 
1957). 

[27] F. Gross, Relativistic quantum mechanics and field theory, (Wiley, New York, 1993). 
[28] J. Collins, Renormalization, (Cambridge Press, Cambridge, 1984). 
[29] G. Baym, Phys. Rev. 117 (1959) 886. 

[30] G.C. Wick, Phys. Rev. 96 (1954) 1124; R.E. Cutkosky, ibid. 1135. 
[31] M. Mangin-Brinet and J. Carbonell, Phys. Lett. B 474 (2000) 237. 
[32] C.R. Ji, Phys. Lett B 322 (1994) 389. 

[33] S. Ahlig and R. Alkofcr, Ann. Phys. (N.Y.) 275 (1999) 113. 

[34] R. Rosenfelder and A.W. Schreibcr, Phys. Rev. D 53 (1996) 3337 & 3354. 

[35] S. Glazek et al, Phys. Rev. D 47 (1993) 1599. 



10 



